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ABSTRACT 


Chemical  lasers  pumped  by  the  reaction  of  atomic  fluorine  with  molecular 
hydrogen  emit  power  from  rotation-vibration  transitions  of  excited  HF  with 
upper  levels  as  high  as  v  =  3.  Collisional  processes  compete  vith  stimulated 
emiesion  for  the  energy  of  the  excited  riF.  A  simplified  analysis  is  presented, 
here  for  intensity,  energy,  and  chemical  efficiency  of  a  class  of  c»ach  lasers. 
Results  are  obtained  in  closed  form.  A  comparison  with  more  exact  computer 
solutions  establishes  the  validity  of  the  analysis  despite  its  simplifications. 

A  comprehensive  parametric  study  examines  the  relative  importance  of  initial 
conditions,  optical  cavity  parameters,  rtnd  rate  coefficients  for  pumping  and 
deactivation  reactions. 
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1.  INTRODUCTION 


This  paper  presents  a  closed-form  analytical  description  of  an  HF 
chemical  laser  pumped  by  the  reaction  of  atomic  fluorine  with  molecular 
hydrogen.  Simultaneous  effects  of  laser  power  extraction,  chemical  pumping, 
and  chemical  deactivation  are  considered  within  the  idealization  of  volume 
uniformity  of  these  processes.  Hence,  the  results  are  most  applicable  to 
lasers  having  premixed  reactants,  as  in  the  pulsed  laser  analysis  of  Aircy.* 
The  analysis  is  also  applicable  to  a  cw  laser  when  the  reactants  are  rapidly 
mixed  and  the  cavity  is  a  Fabry- Perot  resonator.  Modifications  to  the  pulse 

analysis  for  description  of  cw  operation  are  given  in  an  apDendix. 

2 

One  of  the  authors  has  developed  a  comprehensive  computer  program 
that  calculates  the  performance  of  a  chemical  laser  by  numerically  integrating 
a  system  of  nonlinear  equations.  To  complement  such  computations,  we  have 
undertaken  the  present  investigation  with  the  objective  of  deriving  approximate 
closed-form  expressions  for  intensity,  energy,  and  chemical  efficiency. 
Features  of  the  present  simplified  model  and  the  computer  model  arc  sum¬ 
marized  in  Table  I  and  discussed  in  the  following  paragraphs. 

Section  II  presents  the  formulation  of  an  analytical  model  for  a  multi¬ 
level  laser.  We  assume  that,  from  lasing  threshold  to  cutoff  there  is  simul¬ 
taneous  lasing  between  all  pairs  of  adjacent  vibrational  levels.  This  assump¬ 
tion  allows  the  use  of  gain  expressions  to  establish  the  relative  populations  of 
different  Vibration -l  levels.  For  tractability ,  we  restrict  our  attention  to  a 
Doppler -broadened  gain  profile. 

-  1  - 


Table  I.  Comparison  of  Theoretical  Models 


Section  III  presents  the  simplified  chemical  kinetic  model  that  is 

3  4 

employed  in  the  analysis.  Guided  b  «•  results  of  computer  calculations  ' 
employing  a  more  comprehensive  kinetic  model,  we  select  predominant 
reactions  and  treat  most  of  these  as  unidirectional  in  their  dominant  direction. 
For  tractability,  a  constant  temperature  is  assumed.  Consequently,  the 
analysis  is  most  useful  for  mixtures  containing  large  amounts  of  inert  diluent. 

The  multi-level  model  is  specialized  to  the  case  of  an  HF  laser  in  which 
pumping  is  furnished  by  the  reaction 

F  +  H2  -» HF(v)  +  H,  v  =  0,  1,  2,  3 

Expressions  are  derived  that  show  the  effects  of  initial  composition,  tempera¬ 
ture,  and  optical  cavity  output  coupling  on  laser  performance.  Interestingly, 
the  results  are  independent  of  rate  coefficients  for  VV  collisions  between  HF 
molecules . 

Selected  cases  are  compared  with  results  obtained  with  the  more  com¬ 
prehensive  computer  mocW.  The  a^eement  is  seen  to  be  good.  The  effects 
of  the  parameters  are  illustrated  by  consideration  of  numerical  examples. 
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II.  FORMULATION  OF  THE  MULTI-LEVEL  LASER  MODEL 


Rate  equations  for  a  chemically  reacting  system  in  which  lasing  may 
occur  between  adjacent  vibrational  levels  are 

Pllt^  Xch<‘>  +  Xrad«)  -  Xrad(»t 

”  (1) 

•  •  • 

dn(v-) 

p-dt“  =  X^fv,)  -  Xrad(vf  -  1) 

J 

Here  n(v)  is  the  number  of  moles  in  vibrational  level  v  per  unit  'ass,  p  is 

density,  Xch(v)  is  net  production  of  n(v)  by  chemical  reactions,  Xra(j(v)  is 

production  of  n(v)  by  lasing  from  v  +  1  -*v,  and  is  the  highest  vibrational 

level  considered.  In  the  absence  of  lasing  between  a  pair  of  levels,  the 

appropriate  Xracj(v)  *s  zero.  We  assume  that  lasing  is,  in  fact,  occurring 

between  all  the  levels.  All  Xracj's  are  thus  positive  during  a  time  interval 

(tg,  tc),  where  t^  is  a  threshold  time  and  t  is  the  lasing  cutoff  time. 

Justification  for  this  approximation  may  be  fo^nd  in  discussions  of  computer 

3  4 

modeling  of  HF  lasers.  ’  Study  of  cw  and  pulsed  lasers  pumped  by  the 
F  +  reaction  shows  that  the  bulk  of  the  laser  energy  extraction  (>90  %) 
occurs  when  there  is  simultaneous  lasing  on  the  v  =  1 -►  0,  2-*  1,  and 
3  -*  2  bands. 


Preceding  page  blank 
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Equations  (1)  are  easily  inverted  for  the  Xra(j(v) 


W0)  =  Pdtn(0)  -  Xch(0) 


Xrad (1)  =  P^n(°)  +  n(1)]  -txch(0)  +  Xch(D] 


Xrad(vf  -  0  =  pdt[n(0)  +  *  ‘  *  +  n(vf  -  1>1 

'  [Xch<°>+  ‘  *  *  +  Xch(vf  -  1)] 


(2) 


The  last  of  Eqs.  (1)  is  not  needed  in  this  inversion.  Summation  of  Eqs.  (2) 
yields 


V1 


Xrad  = 


xrad(v)  ‘  A1 


-  A- 


v=0 


(3) 


where 


Vf1 


A 


(4a) 


A2=  £(vf  '  v)Xch(v) 

v  =  0 


(4b) 
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The  quantity  Xra(j  *8  central  to  the  analysis,  since  it  is  used  in  the  determination 

of  puL'e  energy,  intensity,  and  efficiency. 

If  we  assume  that  ail  lasing  transitions  have  the  same  wavenumber  to, 

5 

then  the  cumulative  pulse  energy  per  unit  volume  of  gas  E(t)  is  simply 


E(t)  =  hcNAco /  Xrad  dt  =  1 1.  96  o) /Xrad  dt  (5) 

^  t() 

3  - 1 

Here  E  is  in  joule /cm  and  to  is  in  cm 

We  later  use  the  F  +  reaction  for  pumping.  Available  chemical  energy 
from  this  reaction  is  1.40  X  10  J/mole  of  F  atoms.  With  the  notation  (  )q  for 
an  initial  value,  the  chemical  efficiency  n  is,  in  percent, 


h  = 


11.96xl02w  C 

- 5 -  /  *rad  dt 

1.4xlObp(n  )J  rad 

Fo  t0 


(6) 


where  n^,  is  the  mole -mass  ratio  of  atomic  fluorine. 

For  calculation  of  Xrad>  the  n(v)'s  that  occur  during  lasing  may  be  deter 
mined  by  consideration  of  the  gain  required  for  lasing  in  a  specified  optical 
cavity.  For  a  cavity  with  mirror  reflectivities  r^  and  r^,  the  required  gain 
per  unit  length  g  is  controlled  by^ 


r  r  2Lg(v,  J)  a  j 
r 1 r2  1  ’ 


0, 


vf  -  1 


(7) 


where  L  is  the  length  of  active  medium  between  the  mirrors,  and,  for 
simplicity,,  scattering  and  diffraction  losses  are  not  included.  If  rotational 
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equilibrium  at  the  translational  temperature  is  assumed,  the  gain  per  unit 

5 

length  at  line  center  of  a  Doppler -broadened  P-branch  transition  is  given  by 

g(v,  J)  =  grJ  e"^J  "J^6[n(v  +  1)  -  n(v)  e'2J6]  (8) 

Here  the  reference  gain  gr  and  temperature  ratio  6  are  defined  by 

gr  =  2.74  x  1047  <r®rW1/2T'3/2p|^|2,  6  =  ®r/T  (9) 

where  <r  is  a  symmetry  factor,  ®r  is  the  characteristic  rotational  temperature, 
W  is  the  molecular  weight  of.  the  lasing  species,  T  is  the  translational  tempera¬ 
ture,  and  Cni  is  the  matrix  element  of  the  electric  dipole  moment. 

2  5 

The  output  intens  ity  in  W/cm  is  given  by 

I  ~  H.96  «[l  -  (^r2)1/2]xrad/g  (10) 

We  assume  that  in  the  set  of  Eqs.  (8)  the  lower-level  rotational  quantum 

number  J  is  the  same  for  all  lasing  transitions.  This  value  for  J  will  later 

be  chosen  to  maximize  the  efficiency.  We  were  led  to  this  simplification  by 
b  3  4 

previous  analysis  and  numerical  solutions  ’  in  which  most  of  the  power 
(>75  %)  emanates  from  only  two  or  three  transitions  for  each  v  +  1  -*  v  band. 

The  J's  oi  these  dominant  transitions  do  not  vary  much  with  v.  The  use  here 
of  a  single  J  for  alt  levels  and  for  the  entire  pulse  duration  permits  a  simple 
and  reasonably  accurate  total  energy  prediction.  In  consequence,  however,  we 


vH- 


sacrifice  detailed  prediction  of  spectral  content  and  obtain  only  a  single 
dominant  J. 

With  this  dominant  J  assumption,  Eqs.  (8),  with  their  lefthand  sides 
replaced  by  the  g  satisfying  Eq.  (7),  can  be  solved  for  the  n(v)'s.  After 
manipulator ,  the  result  is 


v  =  0,  1,  ■  •  ,  v{ 


IH.  CHEMICAL  KINETICS 


Dominant  kinetic  processes  involved  in  the  class  of  HF  lasers  pumped 
by  the  F  +  reaction  are  now  described.  For  each  reaction  we  define  vibra¬ 
tional-level  distribution  constants  a(v)  such  that 

v, 

^a(v)  =  i  (13) 

For  example,  the  forward  rate  coefficient  k(v)  for  the  pumping  reaction  is 

written  as  k(v)  =  a(v)k(T),  where  k(T)  is  the  overall  forward  rate  coefficient. 

Our  choice  of  reactions,  overall  rate  coefficients,  and  distribution  constants 

is  shown  in  Table  II.  This  choice  is  based  on  analysis  of  computer  solutions 

3  4 

for  flow  and  pulsed  HF  lasers.  ’  Numerical  values  for  the  rates  and  distri¬ 
bution  constants  are  taken  from  Ref.  4. 

Disregarding  backward  rates  for  Reacts.  (R-l)-(R-4)  is  legitimate 
during  lasing  but,  of  course,  is  incorrect  later  when  equilibrium  is  approached. 

Reactions  (R-2)  and  (R-4)  are  the  dominant  VT  processes.  For  React. 

7 

(R-2)  we  use  a  numerical  fit  to  rate  data  of  Airey  and  Fried.  The  distribu¬ 
tion  constants  in  Reacts.  (R-2)  and  (R-4)  are  based  on  a  harmonic  oscillator 
approximation. 

Reactions  (R-3),  (R-5),  and  (R-6)  represent  dominant  VV  processes.  In 
application  of  React.  (R-3)  we  shall  use,  in  the  lefthand  side,  total  concentra¬ 
tion  of  rather  than  H^iO).  The  error  incurred  is  small.  Backward  and 
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forward  reactions  for  HF  -  HF  VV's  are  comparable  in  magnitude.  Therefore, 
we  include  the  forward  reactions  as  React.  (R-5)  and  the  backward  reactions 
as  React.  (R-6).  Nonzero  rate  and  distribution  constants  for  Reacts.  (R-5) 
and  (R-6)  are  given  only  symbolically  in  Table  II.  Later  these  symbols  drop 
out  of  the  analysis.  No  VV  or  VT  processes  that  produce  HF(v)  with  v  >  v^ 
are  included.  These  processes  are  weak. 

A  process  not  included  in  Table  II  is  VT  deactivation  of  HF  by  the 

A 

diluent.  Such  processes  have  rate  coefficients  several  thousand  times 
smaller  than  those  of  React.  (R-2),  and,  unless  the  diluent  concentration  is 
inordinately  large,  they  may  be  disregarded.  As  a  consequence,  the  type  of 
diluent  and  its  concentration,  providing  it  is  sufficient  for  the  isothermal 
approximation  to  apply,  do  not  enter  into  the  analysis. 

On  the  b-sis  of  the  reactions  in  Table  II  we  can  write  Xch(v)  as 


^ch(v)  =  P  |i(v)knj.njj  i  kj_£Fn£jF[a£|F(v  ^  l)n(v  +  1)  -  *^£jp’(v)n(v)  ] 

2* 

+  kFnF[aF(v  +  l)n(v  +  1)  -  ap(v)n(v)] 

+  kH2nH2[aH2*V  +  1*n(v  +  ^  "  at^(v)n(v)] 

+  kfHF[af*v  +  ^n2^v  +  "  2af(v)n2(v)  +  af(v  -  l)n2(v  -  1)J 

+  kbHF^’ab^v  +  +  +  2ajj(^')n(v  +  l)n(v  -  1) 

-  ab(v  -  l)n(v)n(v  -  2)]|  (14) 


where  a(-i)  is  zero,  and  v  -  0,  1,2. 
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In  addition  to  the  reactions,  the  n^  satisfy  atom  conservation  relations 
nH  +  nHF  +  2nH2  =  (®H  +  nHF  +  ZnH^  *  nF  +  nHF  =  *nF  +  aHF  q  <15) 


where  n^  =  n^. 


-H- 


IV.  ANALYSIS 


We  analyze  an  idealized  situation  in  which,  at  t  =  0,  a  homogeneous 
mixture  with  reactant  composition  (n^)  =  0,  (np)  ,  »  an^  (»jjp)  begins 

to  react  in  accord  with  the  kinetic  model  of  Section  III.  Assuming  that  lasing 
threshold  is  reached  immediately,  i.e.,  that  =  0,  we  determine  the  time- 
dependence  of  the  reaction  and  then  calculate  intensity,  shutoff  time,  energy, 
and  efficiency. 

The  composition  of  F  and  H  changes  through  React.  (R-l)  with  rate 
equations 


dtnF  =  “dtnH  =  "kpnFnH- 


We  thus  have 


np  +  rij^  -  (np  n._j) 


Use  of  Eqs.  (15)  and  (17)  in  (16)  gives 


d 

dtnF 


'kpnr[nF  +  (nn2  "  nr)0J 


which  integrates  to 


InTT 


(1  -  k  j )  e 


•  t/r 


1  -  *  j  e 


-t/r 


(16) 


(17) 


(18) 


(19) 
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where 


(3  •  t'‘  b  kpK  -  -f), 


When  =  1,  solution  (19)  is  indeterminate,  and  the  appropriate  form  is 


[1+k<’<V0t]‘1’ 


•i  * 1 


The  time-dependent  composition  can  be  expressed  in  terms  of  a  single 
variable,  the  extent  of  reaction,  defined  by 


From  Eqs.  (15),  (17),  and  (22)  it  can  be  established  that 


nHF  =  (nFJ  J  +  (nHFl 


nF  r  -(nF^  + 


nH  =  ^F^  +  ^nH^0 
nH  =  -,nF>  &  +  ("H,)_ 


and  use  of  these  in  (19)  gives 


1  -  e 

1  -  * je 


We  now  turn  to  a  derivation  of  Xrad-  Use  of  expression  (4a)  for  A  , 
incorporation  of  Eq.  (12)  with  n  t  =  nHF  and  vf  =  3,  and  finally  use  of  Eq.  (23) 
yields 

A  =  m.p(n„)  ^  (25) 

i  i  *  o 

where  the  function  which  is  plotted  in  Fig.  1,  is  given  by 


m 


-  -2J6  ,  -4J6 

3  +  2e  +  e 


1  (> 


+  e-2J6)(l+e-4J6) 


(26) 


Differentiation  of  Eq.  (24)  and  combination  with  Eq.  (25)  gives 


A1  =’mlke  (npn^jjl  -  6)(1  -  «,l) 


(27) 


Substitution  of  Eq.  (14)  into  Eq.  (4b)  gives  A2«  For  brevity  we  en'er 
numerical  values  for  the  a(v)'s  to  obtain 

A2  =  p2jo.  945  knpnH  +  £(kHFnHF  +  kjji^nd)  +  2n(2)  +  3n(3)] 


+  ku  n.,  [0.965  n(l)  +  0.035  n(2)j 
H2  H2 


(28) 


The  HF  —  HF  VV  terms  from  Reacts.  (R-5)  and  (R-6)  have  cancelled  as  a 

result  of  the  assumntions  of  equal  spa  ;g  for  all  vibrational  energy  levels 

3  4 

and  simultaneous  lasing  on  all  bands.  In  computer  modeling,  ’  where 
anharmonic  spacing  is  accounted  for,  these  VV's  have  a  very  small  effect  on 
total  power. 
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Figure  1.  Parameters  m.  vs  J <£'  /T 
o  i  r 


(See  Eqs.  26  and  20. ) 


We  further  simplify  by  setting  0.945  1  and  0.965  n(l)  +  0.035  n(2) 

£*n(l).  Into  this  we  substitute  Eq.  (11)  to  obtain 


2l  nHF  f 

A2  ~  p  jknFnH2  +  m2(kHFnHF  +  kFnF)  +  kH  nH  ] 


g  e(j2-j)6r 

+  J  m4(kHFnHF  +  kFnFJ  +  m5kH^nH, 


B, 


(29) 


'2  2- 


where  the  functions  •  •  •  m^,  which  are  plotted  in  Fig.  1,  are  given  by 


m 


2  =  £(1  +  2e"2J6  +  3e"4J6)(  m3  =  (1  +  e^u)(l  +  e"’JU) 


2J6Wi  .  -4J6 , 


m , 


2J6  -2J6  ,  ... 

e  +e  +  (4/3) 


2J6  -2J6  , 

e  -  e _ -_Z_ 

m.. 


(30) 


Using  Eq.  (23),  we  express  A^  in  terms  of  the  extent  of  reaction 
parameter 


-  kf’  (nFnH,) 

2  n  ( 


(l  -  0(1  -  «iH 


/£ +  *  >  \ 


P'  \  rm 


4  m5oj(l  -  x  jOkH 


*•  * 


m 


(c  4  «  ,)kHp  s-  (1  -  Okpj 


1) 


Fj| 


( 3  j  1 


where 


•P&- 


g  e 
a  =  — 


(J2-J)6 


(nF) 


(32) 


and  ratios  of  rate  coefficients  to  that  for  the  pumping  reaction  are 


--  /ku  /kl 


h„  rH2  j 


kp  =  (kp/k), 


kHF  =  'kHF/k) 


(33) 


These  ratios  of  rates  are  shown  as  functions  of  temperature  in  Fig.  2. 
Substitution  of  (27)  and  (31)  into  (3)  yields  the  central  result 


Xrad  =  kP2|nFnH  j  j^ml  "  "  * 


1^)  -  (■ 


£  ^  *  2  \  _ 

"rn^  +  m5 7(1  '  *l^)kH2 


"m- 

■*lLs;(6  +,2,tm4^ 


($  +  «2)kHF  +  (I  -  eikF 


(34) 


With  J  and  a  set  of  initial  conditions  specified,  quantities  on  the  righthand 
side  are  constant  except  the  timelike  parameter  £.  The  deactivation  terms, 
i.e.,  those  involving  ,  k^p,  and  kp,  all  diminish  Xracj*  We  define  as 
that  value  of  f,  for  which  is  reduced  to  zero.  Lasing  cuts  off  at  %  -  £c> 


and  the  corresponding  time,  obtained  Irom  Eq.  (24),  is  tc> 


The  time  integral  of  X. 


gives  E  and  q,  Eqs.  (5)  and  (6).  It  is 


/ 

yt„=o 


X  ,  dt  =  : — r- 
rad  kp/n^ 


*  xrad(6)  « 


SJVo 
c  0 


pfapMOrij  -  i)|  +  -^—[4  +  (1  +  «2  +m3m5a)  <n(l  -  4)] 


kFf  ,  .  +*l«2>m2  +,1 

^  H6  + - 51 - 

W  ,  ."l*1  +V«!  + 
—  [-m25  + - r 


m^m^a 


<n(l  -  * 


.  kHF  r  ,.«,<•+ Vl(,+,Z)mZ +  tn3nV'l.  ,, 

+  -57  [_-m2^ + - — - <nd  -  ei 


(1  +  *1*2)[(1  +  *  1*  2^m2  +  k 


^(1  -  *  ) 


t n(l  -  *  t4  ) 


As  mentioned  previously,  J  is  chosen  as  that  integer  which  maximizes  t). 

Simple  scaling  laws  are  superficially  indicated  by  Eqs.  (34 ) and  (3 5)  if  one 

neglects,  temporarily,  the  pressure-dependence  of  a.  Output  power  is  pro- 
2 

portional  to  kp  |npn^  \  and  hence  appears  to  scale  as  the  pressure  squared. 

2'q 

Pulse  energy  is  proportional  to  p(np)  and  appears  to  scale  linearly  with  the 
pressure.  Efficiency  appears  to  be  independent  of  the  pressure.  It  will  be 
seen  later  that  the  pressure-dependence  of  a  alters  this  simple  scaling  sub¬ 
stantially. 

For  limiting  cases,  Eqs.  (34)  and  (35)  are  simplified,  and  examination 


of  certain  features  of  the  theory  is  facilitated. 


A. 


THRESHOLD  GAIN  VANISHING 


This  is  the  case  g  =  0,  which  is  approached  as  L  ■♦oor  as  r ^r,  -*  1. 

In  Eqs.  (34)  and  (35)  we  obtain  this  limit  by  setting  a  =  0.  Inspection  of 
Eq.  (34)  shows  that  deactivation  is  less  for  a  =  0  than  for  finite  a.  Hence, 
as  expected,  best  performance  is  predicted  for  the  most  effective  cavity,  i.e., 
for  g  =  0.  With  a  =  0,  maximum  efficiency  occurs  when  J  -*»,  in  accord  with 
the  analysis  in  Ref.  5.  Later  we  shall  find  that  whenever  the  a  terms  are 
nonzero  maximum,  efficiency  occurs  at  «  finite  J. 


B. 


NO  DEACTIVATION 


In  this  limit  k„ 
H2 


=  0  and  the  theory  predicts 


Xrad  =  kPZ(nFnH  )  (mi  '  1)(1  '  ^)(i  '  K  li]  (36) 

'  2  0 

!i /*  ^  for  iij  x  i 

(37) 

1  for  *  j  <  1 


E(t)  =  11.96  up (n 

„  _  11.96  x  iq2 
9  5 

1 . 4  X  103 


F)  (m  -  1)| 
F  0  1 


u(mj  -  l)^c 


(38) 

(39) 


Equation  (37)  predicts  cutoff  when  n^,  or  reaches  zero  and  no  further 
HF(v)  can  be  formed.  The  intensity,  which  is  proportional  to  Xracj,  decreases 
parabolically  with  £  from  its  maximum  value  at  £  =  0.  Maximum  efficiency  is 
achieved  when  nrij  has  its  maximum  value  of  3,  and  this  occurs  for  J  —  co. 


Curiously,  in  this  limiting  case  there  is  no  dependence  on  g.  This  is 
caused  by  idealizations  made  in  the  development  of  Eq.  (12)  and  also  by 
neglect  of  the  time  interval  of  pumping  needed  to  achieve  threshold. 


C.  HYDROGEN-RICH  MIXTURE 


by 


The  case  of  a  large  excess  of  hydrogen  with  fixed  fluorine  is  described 

-*0, '  (n_,)  fixed.  In  this  limit 
F  0 


rad 


,  2,  .2 

ke  (df)o 


(mt  -  1 )( 1  -  4) 


/!+'  2 

-  rsr~ +  m5"ikH- 


(40) 


Ti  =  8.  55  X  10"3  wjfnij  -  1)4  +-“■[4  +  (1  +  *2  +  n^n^o-)  t n(l  -  4)]|  (4  ' ) 


Fast  pumping  and  deactivation  by  predominate.  Because  is 

2 

small  in  comparison  with  kp  and  kj^p  (see  Fig.  2),  a  very  small  *  ^  is  required 

for  this  limit  to  apply. 

D.  FLUORINE-LEAN  MIXTURE 

For  («_,)  —0,  negative  X  .'s  are  predicted.  This,  of  course,  is 

J?  q  rad 

physically  incorrect.  In  reality,  when  fluorine  is  very  scarce,  threshold,  if 
achieved,  is  not  achieved  quickly.  Thus,  the  theory  is  not  suited  to  describing 
this  limiting  case,  because  the  assumption  tg  =  0  is  inappropriate. 

The  preceding  limiting  cases  in  part  illustrate  expected  behavior  and  in 
part  serve  to  point  out  regimes  in  which  underlying  assumptions  of  the  theory 
are  not  fulfilled.  On  the  basis  of  arguments  similar  to  those  for  the  fluorine- 
lean  situation,  we  have  found  that  the  initial  HF  parameter  * ^  should  be 
restricted  to  values  small  in  comparison  with  unity  for  the  theory  to  hold. 
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V.  COMPARISON  WITH  COMPUTER  MODEL 


Development  of  this  closed-form  model  for  an  F  +  laser  was  aided 
by  simplifications  suggested  by  study  of  computer  modeling  of  the  same 
laser.  In  this  section,  we  compare  predictions  of  our  closed-form  model 
with  those  of  the  more  elaborate  computer  model.  Basic  features  of  each 
model  were  summarized  earlier  in  Table  I. 

We  consider  cases  in  which  the  initial  gas  mixture  has  a  composition 
with  mole  ratio  SF^t^tFrHF  =  50: i :  and  a  pressure  of  20  Torr.  A 

large  amount  of  inert,  high-heat-capacity  diluent,  SF^,  ensures  a  nearly 
isothermal  gas  condition.  At  20  Torr,  the  transitions  are  largely  Doppler- 
broadened.  The  length  L  is  fixed  at  10  cm,  and  the  reflectivity  r ^  is  taken 
as  1.  0.  In  anticipation  of  a  parametric  study  over  ranges  of  values  for  T, 

*1*  *2’  ant*  rl  *°  Presented  later,  we  choose  cases  that  provide  a  sampling 
of  these  ranges  for  comparison  with  the  computer  model. 

Values  for  parameters  in  g^,  Eq.  (9),  and  for  o  appropriate  for  HF 
5 

are: 

<r  =  1,  ®r  =  30.  1 6* K,  W  -  20g/mol, 

I  9W  =2.8x10  ^  erg /cm 3,  w=  3790  cm  * . 

Illustrations  of  the  choice  of  a  dominant  J  for  the  present  theory  are 
considered  before  we  go  on  to  the  comparison  with  the  computer  model.  To 
obtain  a  dominant  J,  we  try  several  integer  J  values  until  a  J  giving  maximum 
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il  it  found.  Results  of  such  a  procedure  are  shown  in  Fig.  3.  For  each 
example  the  choice  J  =  7  gives  maximum  efficiency.  Neighboring  J'a  give 
moderately  smaller  efficiencies.  Dominant  J's  were  chosen  in  this  way  in 
all  calculations. 

Table  III  presents  a  comparison  between  the  theory  and  computer  mod¬ 
eling  for  four  cases.  In  the  table,  is  the  time  required  for  all  three 
v  +  1  v  transitions  to  reach  threshold,  Tc  is  the  temperature  at  cutoff  of 
lasing,  and  "output  J"  refers  to  values  of  J  for  transitions  yielding  the  bulk 
of  the  output  energy  for  all  vibrational  bands.  The  results  for  t^  and  Tc 
justify  our  assumptions  of  t^  =  0  and  isothermal  conditions,  respectively. 

The  theory  predicts  values  somewhat  too  large  for  the  dominant  J  and  the 
efficiency.  However,  the  theory  is  accurate  to  within  a  factor  of  i.  3  for  the 
efficiency  for  all  cases. 

On  the  basis  of  the  results  of  Table  III,  the  use  in  the  theory  of  a  single 
cutoff  time  appears  to  be  questionable.  Figure  4,  which  shows  intensity  vs 
time  for  Case  A,  demonstrates  that  the  numerical  calculations  typically  have 
a  long,  low-energy  tail  that  can  be  disregarded.  Thus,  the  large  variation 
in  cutoff  times  is  not  as  serious  as  first  appears.  The  overall  agreement  of 
pulse  shape  shown  in  Fig.  4  is  quite  good. 


Table  III.  Predictions  of  Theory  Compared  with  Those  of  Computer  Model 


Case 

Specified  Parameters 

B 

C 

D 

T,  °K 

200 

300 

400 

300 

*1 

0.  5 

1 

2 

1 

*2 

0 

0 

0 

0.  1 

rl 

0.8 

0.  6 

0.4 

0,  6 

Computer  Model  Results 

tQI  (isec 

0.  12 

0.  14 

0.20 

0.  54 

(1  0)  tc,  psec 

21.9 

12.  9 

6.  1 

11.4 

(2  —  1)  tc>  psec 

36.  1 

18.  5 

7.2 

16.8 

(3  —  2)  t  ,  psec 

19.5 

6.4 

2.4 

5.9 

Outp  it  J 

4,  5 

4,  5,  6 

5,  6 

5.6 

T  ,  "K 
c’ 

210 

311 

410 

311 

r\ 

35.4 

27.  8 

13.8 

24,7 

Results  of  Theory 

t  ,  psec 

24.  1 

12.  9 

4.  9 

11.  5 

Dominant  J 

5 

6 

7 

6 

9 

36.  4 

32.  1 

18.  1 

29.  5 

Figure  4.  Output  Intensity  vs  Time  for  Case  A  of  Table  TIT 
(Assumptions  for  the  computer  model  are  listed  in  Table  I.) 


VI.  PARAMETRIC  STUDY 


The  theory  is  used  to  study  the  dependence  of  laser  performance  on  the 
parameters  —  pressure  p,  temperature  T,  initial  fluorine -to -hydroden  ratio 
initial  HF-to -fluorine  ratio  k^,  and  mirror  reflectivity  r j  for  fixed  values 
of  length  L  =  10  cm  and  reflectivity  =  1.0.  The  relative  importance  of 
pumping  and  deactivation  reaction  rates  is  assessed  with  the  aid  of  reaction 
rate  parameters  defined  as  follows: 


d 

P 


=  (m 


1 


-  1)£ 


=  (1  +  * 
m3 


d„  =  —  <m0£  +  * 


m 


1 


2  +  m3m5o')  ln(l  -  £)] 

[  (1  +  *2^m2  +  K^m^m^ft]  fn(l  - 


(42) 


(43) 


(44) 


For  the  form  depends  on  the  value  of  «j.  For  4  1 


dHF^l  *  "  m,  rm2^  +  '  *1^  V1  +  X2) 


[(1  +  «2)m2  +  m3m4tt]  in0  -  £)  -  («j  -  *  j)  (1  +  KiKz) 


x  [  ( 1  +  *  i*Z)m2  1  K  jm3m4Q3  -  *'jC)| 


(45a) 
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and  for  *j  =  1 


kj|p  |  '  > 

djjp(«i  =  1)  =  j" m2^  "  ^  *  *2^m2^  *  K  2}  +  m3m4a^  j  _  £  (45b) 

-  [2(l  -  *2)m2  +  m3m4“J  *n(l  "  £) 

In  terms  of  the  d's,  pulse  energy  and  efficiency  are  given  by 

E(t)  =  11.96  »P(»F)0(<ip  +  d^  +  dF  +  dHf,)  (46) 

T1  =  8.  55  X  10'3  «|dp  +  dH  +dF  +  dHFJ  (47) 

Thus,  the  d's  allow  us  to  assess  the  importance  of  the  deactivation  mechanisms 

relative  to  each  other  and  to  the  pumping  mechanism  dp. 

Results  of  a  parametric  study  for  initial  gas  composition  with  mole  ratio 

SF^:Hj,:F-HF  =  50:  If  K1K2  anc*  P  =  Torr  are  presented  in  Table  IV  for 

combinations  of  T  =  200°K,  300°K,  400°K  and  r^  =  0.  8,  0.  6,  0.;4  for  the  cases 

=  0.5,  =  0;  =  1.0,  *2  =  *1  =  2.:  0,  ^  -  0..  and  =  l,  Q,  -  0.1. 

Table  IV  shows  that  is  negligible  in  comparison  with  dp  and  d^p. 

This  is  consistent  with  the  small  values  of  kTI  i  comparison  with  those  of 

H2 

kp  and  k^p  in  the  temperature  range  of  interest  (see  Fig.,  2).  Considering  the 
main  deactivators  we  see  that,  with  few  exceptions,  d^p  is  larger  than  dp.: 
However,  when  =  2  and  T  =  300  and  400°K,;  excess  tluorine  and  high 
temperature  cause  dp  to  predominate  Throughout  the  table,  increasing 


Table  IV.  Parametric  Study 


T.’K 

'i 

J 

H»*e 

d. 

■s 

"r 

ilcm1 

*1 

■j  ■  0.  9 

*2*0 

209 

0.1 

S 

0.005 

24  15* 

1.300 

tl 

f.  104** 

2.555° 

2. 170' 1 

1.  1*1° 

54.55 

0.4 

5 

0.  777 

11.75 

1.340 

1.051° 

5.416° 

2.  Oil'1 

7.651° 

34.04 

0.4 

5 

0. 741 

19.  14 

1.270 

5.047° 

4. 175° 

2.577° 

7.612° 

31.2* 

MO 

O.t 

7 

0.017 

15.  64 

1.550 

2.752° 

5. 555"* 

1.114** 

5.421° 

42.02 

3.C 

4 

o.on 

15.35 

1.443 

2.000° 

1.  745° 

1.720*1 

5.6*5° 

40.44 

0.4 

4 

0.007 

13.51 

1.424 

2.545° 

5.  155° 

1.554° 

5.665° 

30.  OS 

400 

o.i 

• 

0.  047 

13.41 

1.523 

4.5*5° 

5.074° 

i.  sir1 

5.645° 

45.22 

0.4 

7 

C.  045 

13.03 

1.442 

4  »0»° 

4.5*5° 

1.204“* 

4.  575° 

40. 74 

0.4 

7 

0.024 

11.27 

1.410 

4.2*5° 

7.  5*5° 

1.53l“* 

4.261° 

30.  14 

V 

■**• 

209 

0.1 

4 

0. 422 

17.70 

1.  124 

4.457° 

4.  764° 

2. 104*  * 

1.265° 

27.74 

0.4 

3 

0.407 

14. 44 

1.047 

4.515° 

4  >15° 

1  005* * 

MU*1 

24.04 

0.4 

S 

0.  502 

15.01 

1  004 

4  455° 

5.  161° 

2.  040*  * 

1.667° 

24.23 

MO 

0.0 

7 

0,747 

13.90 

1  204 

1  155° 

5.  545° 

t.oor* 

1  055° 

34.52 

0  4 

4 

0.755 

12.91 

1  209 

1  175° 

5.550° 

1.555° 

5.  556° 

32.00 

0.4 

4 

0.  720 

tl.40 

1.  171 

1.652° 

7.  152° 

1.442’* 

5  626° 

M.  II 

400 

0.1 

1 

0  015 

12.43 

l  M7 

2.515° 

1  >55° 

1.  32 1“* 

7.515° 

55.2* 

0.4 

7 

0  004 

11.04 

1  227 

1.555° 

•  154° 

1.210** 

7  567° 

32.02 

0  4 

7 

0  774 

9,  99 

1.  104 

..755° 

1.055° 

1. 334*  * 

6  «55° 

M.  55 

.,.10 

«j  •  0 

200 

0  I 

4 

0  300 

11  79 

0  721 

1.  156° 

4  *66° 

1.001’* 

1.545° 

10.  S3 

0  4 

4 

0.  MO 

10  17 

0.  470 

1.567° 

4  566° 

1  000’* 

1  525° 

17  03 

0.  4 

5 

0  377 

10.04 

4  504* 1 

1  451° 

4  670° 

1  040* * 

1.  444° 

14.  00 

MO 

0  • 

7 

0  454 

7  340 

7  4T5*1 

2.525° 

5.552° 

5.570° 

1.2*1° 

21.07 

0  4 

6 

0  440 

7.042 

7.  192’ 1 

2.554° 

5  61*° 

5  602° 

1  172° 

>5  55 

0  4 

4 

0  439 

4  440 

7  049* 1 

2.511° 

5  555° 

5.275° 

1  115° 

11  44 

400 

0  • 

• 

0  444 

5  700 

7  447* 1 

5  *55° 

7.  150° 

5.501° 

1.247° 

20.41 

0  4 

7 

0  440 

5  490 

7  014" 1 

>  750° 

7.  115° 

5.  525° 

t.  *55° 

10  25 

0  4 

7 

0  449 

4  0«0 

4  854' 1 

1.  *75° 

1.727° 

5.  555° 

•  115° 

10.  00 

temperature  causes  djjp  to  decrease  in  magnitude  like  the  E^p  curve  in 
Fig.  2. 

Increasing  r.  causes  moderate  increases  in  tj,  £  ,  and  t  .  As  r. 
increases  the  dominant  7  increases,  as  shown  in  Ref.  5.  Here,  however, 
the  increase  in  J  is  never  greater  than  one  when  r^  varies  between  0.  4  and  0.  8. 

The  effect  of  initial  HF  is  shown  by  a  comparison  of  the  group  of  cases 
for  =  1.0,  *2  =  0.  1  with  corresponding  cases  for  * ^  =  1.  0,  * ^  =  0.  In  each 
case,  initial  HF  causes  a  drop  in  efficiency;  the  effect  is  greatest  for  T  =  200 °K. 

For  most  cases  in  Table  IV  the  efficiency  is  largest  when  T  =  400  °K  and 
smallest  when  T  =  200°K.  The  differences,  however,  are  less  than  might  be 
anticipated  from  the  kj^p  and  kp  curves  in  Fig.  2.  Indeed,  one  might  conclude, 
reasoning  solely  from  those  curves,  that  lasing  could  not  occur  at  T  =  200 °K. 
However,  even  at  200°K  and  with  initial  HF  present,  i.e.  ,  *2  -  0.  1,  the 
deactivation  processes  are  stronger  but  are  not  dominant.  This  conclusion 
stems  from  a  comparison  of  jd^  +  dp  +  dj^p  |  with  dp. 

An  exceotion  to  the  main  trend  with  temperature  is  the  group  of  cases  for 
-  2.  0,  *2  =  0,  where  the  largest  efficiencies  are  found  at  300°K.  The  cause 
of  this  behavior  is  the  change  in  tempe  rature -dependence  of  the  deactivation  for 
this  fluorine-rich  case,  m  which  dp  and  dpjp  are  of  nearl,  equal  importance. 

These  trends  of  efficiency  vs  temperature  prompted  us  to  investigate  an 
extended  range  of  temperature  for  the  case  *j  =  1,  k  ^  -  0,  =  0.  8.  For 

this  case  it  was  found  that  efficiency  peaks  at  400°K  and  decreases  monotoni- 
caily  for  higher  temperatures. 
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Throughout  the  study  we  have  found  that  the  dominant  J  corresponds  to 
the  abscissa  interval  in  Fig.  1  for  which  ra^  begins  to  rise  steeply.  The 
interval  is  indicated  by  vertical  dashed  lines  in  Fig.  1.  Low  temperatures 
favor  the  right  side  of  the  interval,  high  temperatures  the  left.  Thus,  larger 
deactivation  at  low  temperature  is,  in  part,  countered  by  larger  m^,  as  shown 
by  Eqs.  (43)-(45)  for  the  d's.  These  countervailing  effects  result  in  the 
existence  of  an  optimum  temperature  for  maximum  efficiency.  This  m^  effect 
is  related  to  the  temperature -dependence  of  the  gain  that  makes  molecular 
lasers  most  efficient  at  low  rotational  temperatures  when  deactivation  processes 
are  not  considered. 

The  effect  of  an  increase  in  the  initial  fluorine -to-hydrogen  ratio  is 

shown  by  changes  from  group  to  group  in  the  first  three  groups  of  data  in 
Table  IV.  Efficiency  is  somewhat  larger  when  is  in  excess,  as  shown  by 
the  =  0.  5  results.  Decreased  efficiency  when  =  2  is  caused  by  our 
fluorine-centered  definition  of  efficiency,  Eq.  (6).  We  can  obtain  another  view 
by  examining  trends  of  pulse  energy  and  pulse  duration  with  For  the  case 
T  =  300°K,  *2  =  0»  rir2  =  8,  we  have  carried  out  calculations  to  supplement 

those  of  Table  IV.  Pulse  energy  increases  with  «tj,  rapidly  at  first,  when 
there  is  excess  and  then  more  slowly  for  greater  than  stoichiometric, 
where  =  1.  Pulse  duration,  after  a  brief  plateau  for  small  decreases 
monotonically.  Shortening  of  the  pulse  with  increasing  xj  is  to  be  expected 
because  of  the  increased  pumping  rates. 

Until  now  the  examples  have  had  a  pressure  of  20  Torr,  The  effect  of 
changing  the  pressure  is  of  interest,  because  at  low  enough  pressures  lasing 
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VII.  CONCLUSIONS 


A  theory  providing  a  closed -form  solution  to  rate  equations  for  a 
chemically  pumped,  multilevel,  molecular  laser  has  been  presented.  Appli¬ 
cation  is  made  to  a  pulsed  HF  laser  pumped  by  the  reaction  F  +  H2~*HF(v)  +  H. 
Application  to  a  cw  laser  is  discussed  in  an  appendix. 

The  analysis  predicts  laser  intensity  vs  time,  pulse  energy,  pulse 
duration,  and  efficiency.  Because  prediction  of  particular  lasing  transitions 
is  avoided,  spectral  information  is  limited  to  prediction  of  a  dominant  J. 
Relative  importance  of  pumping  and  deactivation  rates  is  identified.  To  the 
approximation  of  a  harmonic  spacing  of  vibrational  energy  levels,  it  is 
found  that  laser  performance  predictions  are  independent  of  rate  coefficients 
for  HF  -  HF  VV  transfer.  Good  agreement  is  found  in  comparison  of 
theoretical  predictions  with  those  of  a  comprehensive  computer  model  for 
cases  chosen  to  explore  a  wide  range  of  parameters  but  restricted  to  avoid 
gross  violations  of  the  theoretical  assumptions.  A  parametric  study  is 
presented  that  demonstrates  the  capability  of  the  theory  for  economical 
prediction  of  a  large  number  of  trends. 

Extension  of  the  theoretical  approach  to  other  chemical  systems  seems 
desirable  and  feasible.  The  authors  have  obtained  preliminary  results  for 
chain-reaction-pumped  lasers.  Extension  of  the  theory  for  more  detailed 
predictions  of  laser  performance  also  appears  possible  but  only  at  the 
expense  of  added  complexity,  which  may  prove  to  be  uneconomical  in  com¬ 
parison  with  computer  modeling. 
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APPENDIX.  APPLICATION  TO  CW  LASERS 


We  consider  a  steady  isothermal  flow  in  the  x  direction  with  the  laser 

cavity  axis  transverse  to  the  flow,  and  we  assume  a  Fabry-Perot  cavity. 

With  a  flow  of  constant  velocity  u,  the  independent  variable  t  of  the  pulse 

theory  can  be  replaced  with  x/u.  The  equivalent  to  cumulative  pulse  energy, 

5 

Eq.  (5),  in  a  flow  system  is  the  cumulative  output  power,  in  watts, 

t 

P(t)  =  11.  96  u(Au)  J  Xrad  dt 
l0 

where  A  is  the  flow  cross-sectional  area  and  the  integral  is  given  by  Eq.  (35). 
Efficiency,  Eq.  (6),  and  intensity,  Eq.  (10),  are  unaltered. 
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